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We study the zero-temperature grand potential of a three-dimensional superfluid made of ultracold 
fermionic alkali-metal atoms in the BCS-BEC crossover. In particular, we analyze the zero-point 
energy of both fermionic single-particle excitations and bosonic collective excitations. The bosonic 
elementary excitations, which are crucial to obtain a reliable equation of state in the BEC regime, 
are obtained with a low-momentum expansion up to the forth order of the quadratic (Gaussian) 
action of the fluctuating pairing held. By performing a cutoff regularization and renormalization of 
Gaussian fluctuations, we find that the scattering length ob of composite bosons, bound states of 
fermionic pairs, is given by as = (2/3) as, where as is the scattering length of fermions. 

PACS numbers: 03.75.Ss 03.70.+k 05.70.Fh 03.65.Yz 


I. INTRODUCTION 


Some years ago the crossover from the weakly-paired 
Bardeen-Cooper-Schrieffcr (BCS) state to the Bose- 
Einstein condensate (BEC) of molecular dimers was ex¬ 
perimentally achieved using ultracold fermionic alkali- 
metal atoms in a three-dimensional (3D) configuration 
M- These remarkable experimental results have been 
investigated by several theoretical groups (for a recent 
review see i)- In a recent experiment the BCS-BEC 
crossover has been analyzed also in a quasi-2D atomic gas 
i, where the experimental data are not in full agreement 
with theoretical predictions 00 - In a very recent paper 
[§) we have studied Gaussian fluctuations of the zero- 
temperature attractive Fermi gas in the 2D BCS-BEC 
crossover showing that they are crucial to get a reliable 
equation of state in the BEC regime of composite bosons, 
bound states of fermionic pairs. In particular, per¬ 
forming dimensional regularization and renormalization- 
group analysis, we have found the remarkable analytical 
formula as = af/(2 1 / 2 e 1 ' /4 ) which connects the scat¬ 
tering length as of composite bosons to the scattering 
length of of fermions 0. 

In this paper we investigate the relationship between 
as and of in the 3D BCS-BEC crossover. We study 
the low-momentum expansion, up to the fourth order, 
of the quadratic action of the fluctuating pairing held 
finding that it gives an ultraviolet divergent contribution 
of the Gaussian fluctuations to the grand potential. We 
perform a cutoff regularization of this zero-point energy 
in the BEC regime and obtain the simple but meaning¬ 
ful analytical formula as = (2/3) of- This result is in 
good agreement with other beyond-mean-held theoretical 
predictions: as — 0.75 <zf based on a diagrammatic ap¬ 
proach foj, as — 0.60 clf derived from a four-body anal¬ 
ysis M and also from Monte Carlo simulations [TlJ , and 


as — 0.55 of obtained with convergence factors 00. 
We stress, however, that our result is fully analytical, 
contrary to all other beyond-mean-held predictions 0- 
[l3l | , and it is based on a transparent cutoff regularization 
and subsequent renormalization of bare physical param¬ 
eters. 


II. THE MODEL 

We consider a three-dimensional Fermi gas of ultra¬ 
cold and dilute two-spin-component neutral atoms. The 
atomic fermions are described in the path integral formal¬ 
ism by the complex Grassmann helds ip s {r,T), ip s (r,r) 
with spin s = (t,^) 0, 0. The partition function 
Z of the uniform system at temperature T, in a three- 
dimensional box of volume V , and with chemical poten¬ 
tial fi can be written as 

Z = J DbPsi'&s] ex p|“^ ‘S’l > (!) 


where 


Chf3 


s = 


dr I d 3 r Jgf 

Jv 


( 2 ) 


is the Euclidean action functional and «£? is the Euclidean 
Lagrangian density, given by 


J? = 1ps 


h 2 

hd T -V 2 — n 

2m 


ips + gip ( 3 ) 


where g is the strength of the s-wave inter-atomic cou¬ 
pling (g < 0 in the BCS regime) 0, 0. Summation 
over the repeated index s in the Lagrangian is meant 
and /3 = l/(ksT) with ks Boltzmann’s constant. 
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The Lagrangian density 2zf, which is quartic in the 
fermionic fields, can be rewritten as a quadratic form by 
introducing the auxiliary complex scalar field A(r,r) so 
that El 111: 


Z = J V[ip s gip s ]V[A,A] exp— 
where 


S e (ip s ,i>s, A, A) 


(4) 


SeW’ s ,4>s, A, A) = [ dr [ d 3 r ZZ e (i/ s ,ip s , A, A) (5) 
Jo Jv 

and the effective Euclidean Lagrangian density 
22? e (V>s)Vis! A, A) reads 


^e=$s 


h 2 

hd T -V 2 — g 

2 m 


i> s +A V’i V’t+AV’t ^-i - 


]AJ! 

9 

( 6 ) 

In this paper we investigate the effect of fluctuations of 
the gap field A(r, t ) around its mean-field value Ao which 
may be taken to be real. For this reason we set 


A(r,r) = A 0 + r?(r,r) , 


(7) 


where 7?(r, r) is the complex paring field of bosonic fuc- 
tuations fl4Ml(l. 


III. MEAN-FIELD 

By setting g(r,t) = 0, and integrating over the 
fermionic fields 7/> s (r, t) and , 0 s (r,f) one obtains the 
mean-field (saddle-point and fermionic single-particle) 
partition function (if! ilGI 

Z m f = exp = exp{-(3n m f} , (8) 


where 

S m f 


h = -Tr [In (Go 1 )]- PV 


A 2 
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= - ^ [2 In (2 cosh (/3E sp (k)/ 2)) - /3(e k - g)} 

k 

A 2 

- , 

9 


(9) 


with efc = h 2 k 2 /(2m), 

G - 1 = ( Mr~£;V 2 -H 


Ao 

H 2 V72 


Aq hd T + + /x 


( 10 ) 


the inverse mean-field Green function, and 

E sp (k) = yj (e fc - g) 2 + A 2 (11) 

the energy of the fermionic single-particle elementary ex¬ 
citations. 


At zero temperature (T = 0, i.e. /? —► + 00 ) the mean- 
field grand potential f l m f becomes 


fim/ = - X! ( E sp( k ) “ + 9) - V— 

1, ^ 


( 12 ) 


The constant, uniform and real gap parameter Ao is ob¬ 
tained by minimizing with respect to A 0 , namely 


on 


mf 


OAn 


= 0 . 


In this way one gets the familiar gap equation 
1 _ 1 1 


(13) 


(14) 


where the ultraviolet cutoff A is introduced to avoid the 
divergence of the right side of Eq. (THT) in the continuum 
limit —y V J d 3 k/(2n) 3 . One can express the bare 

interaction strength g in terms of the physical s-wave 
scattering length of of fermions through [id . El,® 


47 xh 2 aF 


_ 1 _ 1 V 1 
~~q ~ V ^ 2el ' 

y | k| < A 


After integation over momenta Eq. m reads 


1 


Anti 2 aF 9 2i t 2 H 2 


A . 


(15) 


(16) 


For this expression one deduces that while g is always 
negative ap can change sign. In particular, in the weak- 
coupling BCS limit, where g —> 0~, the first term on the 
right of Eq. m dominates and of = mg/(Anti 2 ) —> 0 . 
In the strong-coupling BEC limit, where g — 00 , the 
second term on the right of Eq. ra dominates and of = 
7r/(2A) —> 0+ when A is sent to infinity [2p| . 

Inserting Eq. (1151) into Eq. m we obtain the regu¬ 
larized gap equation 


m 


= 1 V 
v 


Anti 2 a f 


1 


V ^ \2E s Jk) 

|k|<A V spV ' 


1 

2 T k 


(17) 


where one can safely take the limit A —> + 00 , finding the 
energy gap Ao as a function of the chemical potential 
g and the scattering length of- We stress that in the 
BCS limit, where aF -> 0”, the chemical potential g is 
positive and g/ Ao — > + 00 . At unitarity, where of —>■ 
± 00 , Eq. m gives g/Ao = 0.8604 showing that the 
chemical potential g is still positive. In the BEC regime, 
where «f —> 0 + , the chemical potential becomes negative 
and it is given by 


ti 2 


V = 


2 ma 2 F 


lmo| 2 
+ A ti 2 0 ’ 


(18) 


while g/A q — 00 . Notice that ep = h 2 /(ma 2 F ) is the 
binding energy of the atomic dimers (composite bosons) 
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which are formed at unitarity E 3 , EM E 3 l and clearly 
[i = —es/2 in the deep BEC limit. 

By using Eq. (fl5l) (with A —> +oo) with Eq. (fT^l) we 
obtain the zero-temperature regularized mean-field grand 
potential of fermionic single-particle excitations 

= ~ E (E'pW - + » - 2ir) “ V 4 • 

(19) 

In the BCS limit, where ar —> 0 , the energy gap Ao 
goes to zero and the mean-field grand potential becomes 
that of a non-interacting Fermi gas, namely 

i™ 

In the BEC limit, where clf —> 0 + , the mean-field grand 
potential reads E 3 


the sound velocity. Both A and c s depend on the chemi¬ 
cal potential p, and the energy gap Ao- In particular, one 

finds EF 


A = 


-li£* BCS limit 
1 unitarity 
BEC limit 


(25) 


and 


b 

= ( b 


BCS limit 
unitarity 

BEC limit 

8 |m| 


(26) 


Integrating over the bosonic fields 77 (g) and 77 (g) in Eq. 
, at zero temperature we find the grand potential 

= - lim \ In (Z) = Q m f + f l g (27) 

/3-H-oo ft 


Vmf = -V 


/ 2m\ 


3/2 


A 4 


2567r \ H 2 J \n\ 3 / 2 


( 21 ) 


We point out that in this limit both Ao and |/z| (/r < 0) 
diverge but \fi\ diverges faster and the grand potential 
H m f goes to zero. 


IV. GAUSSIAN QUANTUM FLUCTUATIONS 


where fl m f is given by Eq. m, while fl g reads 

= 9 El E co i(q) . (28) 

" |q|<A 

This is the zero-point energy of bosonic collective exci¬ 
tations El El [ 20 | . Also here the ultraviolet cutoff A 
is introduced to avoid the divergence in the continuum 
limit Eq~>Vf d 3 q/(27r) 3 . 


We now consider the effect of quantum fluctuations, 
i.e. in Eq. 0 we allow rj(r,t) ^ 0. Expanding the 
effective action S e (ip s , ip s , A, A) of Eq. 0 around A 0 up 
to the quadratic (Gaussian) order in r](r,t) and fj(r,t) 
one finds 

Z = Z mf J £>[77,77] expj-E^lMj ) (22) 

where 

Sg(V,V) = ^J2 {v(q),V(~Q)) M (<2) ( fj{-q) ) ( 23 ) 

is the Gaussian action of fluctuations in the reciprocal 
space with q = (q, iv m ) the 4-vector denoting the mo¬ 
menta q and Bose Matsubara frequencies v m = 27rm//3. 
The 2x2 matrix M(g) is the inverse fluctuation propa¬ 
gator, whose non trivial dependence on q can be found in 
Ref. [H, El • The energy E co i (q) of the bosonic collective 
excitations can be extracted from M (q) El, EME3 and 
it is given by 


d^col id) 


yE (A e q + 2 m c 2 ) 


(24) 


V. SCATTERING LENGTH OF COMPOSITE 
BOSONS IN THE BEC LIMIT 


Expanding Eq. (051) in powers of A [20] we find 


Qg 

V 


___A 5 + 3 A 3 

407r 2 mA 1 / 4 127 t 2 A 1 / 2 

to 3 c 4 . 8 m 4 c 3 1 , 

4tt 2 \ 3 / 2 15n 2 H 2 X 2 V 


(29) 


In this equation the first two terms are truly divergent; 
the third term, despite being oc A is indeed convergent, 
as we are going to show shortly and the fourth term is 
subleading in the BEC limit. 

We then want to regularize the A 5 and A 3 terms and 
this can be done by redefining the bare parameters in the 
mean-field grand potential in eq. ED, introducing their 
renormalized counterparts. One then finds, as shown in 
detail in Appendix,, that the following regularization pat¬ 
tern removes the divergencies at the leading order: 



A 0 + 

ImI - 


16V2 IaiJ ( h 2 \ 
5 -Ajr \2m J 

128 IHJ ( ALA 5 

97T 3 Aq y 2 m J 


5 A 5 
A 5 


(30) 


where e q = h 2 q 2 /(2m) is the free-particle energy, A takes 
into account the first correction to the familiar low- 
momentum phonon dispersion E co i(q) ~ c s hq, c s being 


having used the deep-BEC relations in Eqs. ED, ED 
and, as previously discussed, the fact that in the strong- 
coupling BEC limit the cutoff A can be obtained from 
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Eq. (fl6l) and it reads 



(31) 


We now analyze the leading convergent contribution, 
i.e. the A term of Eq. ED; even though it seems diver¬ 
gent, being proportional to A, actually it is not because 
in the BEC limit cf goes to zero faster than 1/A. 

Using Eqs. (fT51) . Eol) . (ESI) . (I5T1) in the deep BEC 
regime we have \/i\ = h 2 /[2ma 2 F ), A = 1/4, me 2 = 
Aq/(8|/z|), and the A term of Eq. (ESI) becomes 


n B = -v 


a 

256tt 


2 mA 3/2 Ag 

H 2 ) H ^ 72 ’ 


(32) 


with a = 2. It is important to stress that Eq. m is 
formally the same formula found by Hu, Liu, and Drum¬ 
mond [12| and also by Diener, Sensarma, and Randeria 
[bit by using a different regularization procedure based 
on convergence factors. They have determined numeri¬ 
cally the parameter a finding respectively a = 2.5 fl2j 
anda = 2.61 [13], while here we derive a = 2 analytically. 

In the BEC limit, taking into account Eqs. m and 
(1521) the total grand potential is given by 


+ fig —U 


(1 +«) 

256?r 


2 mA 3/2 Aq 

w J ’ 


(33) 


with (1 + a) = (1 + 2) = 3. Here the grand potential Q 
depends explicitly on both and A 0 . Consequently, the 
total number N of fermions must be calculated as follows 


l S / 1 /y ! A 0 V 5A 0/y (I 

and the number density n = N/V reads 

= (1 + «) ( 2mA 37 ' A ° 

167T \H 2 ) l/^l 1 / 2 


(34) 


(35) 


Notice that in Eq. (1551) the second term accounts for a 
variation of the order parameter Ao as a function of 
instead of considering it as an independent variable. The 
same procedure was adopted in Refs. mill , while in 
the well-known Nozieres-Schmitt-Rink scheme this term 
is absent. As stressed in Refs, mill! this second term 
cannot be neglected as it becomes more and more rele¬ 
vant in the BEC regime we consider here, becoming of 
the same order of the first term. 

Eq. (1351) shows that, at fixed number density n, in 
the BEC limit, where both \fi\ and Ao go to infinity, one 
has |//| ~ Aq and, from Eq. (1351) it follows that H goes 
to zero. To obtain Eq. (1551) we have used Eq. ([55]) 
but also Eq. ED, which immediately gives dAo/d[i = 
2h 2 / (ma^Ao) ~ 4|/x| / Ao - Taking into account Eq. (1551) . 
the equation (fl8l) for the chemical potential in the BEC 
limit can be rewritten as 


h 2 

2 ma? F 


7 rh 2 a F 

-TV” —\ n > 

to (1 + a) 


(36) 


where the second term is half of the chemical poten¬ 
tial hb = 47 t H 2 aBnB/niB of composite bosons of mass 
ms = 2m, density ub = n/2 , and boson-boson scatter¬ 
ing length 


dB 


2 2 

——— -a F = -a F 

(1 + o) 3 


(37) 


This result is in good agreement with other beyond- 
mean-field theoretical predictions: as — 0.75aF of Pieri 
and Strinati 9|, as — 0.60a^ of Petrov, Salomon and 
Shlyapnikov [Tcij (and also Astrakharchik, Boronat, Ca- 
sulleras, and S. Giorgini 0), and cib — 0.55 a F of Hu, 
Liu and Drummond [121 ] (and also Diener, Sensarma and 
Randeria 0)- On the other hand the mean-field result 
is quite different, namely as =2 a F 0. As stressed in 
the introduction, the beauty of our result, Eq. J57D, is 
that it is obtained fully analytically, cont rary to all other 
beyond-mean-field predictions @, [lfl, El Ef . 


VI. CONCLUSIONS 

Starting from a theory of attractive fermions, perform¬ 
ing cutoff regularization plus renormalization of Gaus¬ 
sian fluctuations, we have obtained a remarkable formula 
between the scattering length as of composite bosons 
and the scattering length aj? of fermions. This formula, 
as = (2/3) of, is based on the two-body scattering the¬ 
ory which gives in the deep BEC regime an explicit re¬ 
lationship between the cutoff parameter and the s-wave 
scattering length. Our approach, limited to the quartic 
term in the low-momentum expansion of bosonic collec¬ 
tive excitations, is fully reliable in the BEC regime but 
it cannot describe the entire 3D BCS-BEC crossover. In 
fact in the BCS region, where the chemical potential /i 
is positive, the sign of the coefficient A of the collective 
spectrum is negative (pair instability) and further terms 
must be included in the momentum expansion. 
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Appendix: Renormalization in the BEC limit 

By using Eqs. E5l) . E6l) and (1511) we rewrite the diver¬ 
gent terms in the Gaussian grand potential in Eq. ED 
as a function of Aq and ^ 


flf) = 

9 384 




2 TO 

H? 


(A.l) 
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fi (A 5 ) = y V^! 

9 640 




(A.2) 


The tree-level mean-field contribution in Eq. ED is a 
function of the bare parameters Ao and |/x|, which we 
now renormalize as follows: 


for a suitable choice of the counter-terms. The renormal¬ 
ized mean-field grand potential now reads, to the leading 
order: 


Ao —* A/{ = Ao + S A 
M —> \b-\r = H + 6/J. 


(A.3) 

I 


nZf = -v- 


256?r 


2 m \ 3/2 

h 2 ) 


4 

Wr 


-V ■ 


256?r 


2 m\ 

a? J 


3/2 


A 4 
A *0 


I# H = 


A 3 

4 Aio S A- 


3 A 4 




■S[i 


(A.4) 


and in order to find the renormalization pattern in Eq. 
(1301) we simply equate the new terms in the renormalized 
grand potential in Eq. (IA.4[) to minus each divergent 
term, i.e.: 


-v — 

2 256tt 


2m\ 3/2 A 4 (A 3) 


(A.5) 


-4V 


1 

256t r 



-^|<5A = -Qf) 


(A.6) 


solving for the counter-terms one obtains 5/j, = 

yzA 0 

and = fA - Eq - EUl) readily follows reinstating A 
by means of Eqs. E5l) . ESI) and (ED- 
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